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ABSTRACT 

We study the kinematics of satellites around isolated galaxies selected from the Sloan 
Digital Sky Survey (SDSS) spectroscopic catalog. Using a model of the phase-space 
density previously measured for the halos of ACDM dark matter cosmological sim- 
ulations, we determine the properties of the halo mass distribution and the orbital 
anisotropy of the satellites as a function of the colour-based morphological type and 
the stellar mass of the central host galaxy. We place constraints on the halo mass and 
the concentration parameter of dark matter and the satellite number density profiles. 
We obtain a concentration-mass relation for galactic dark matter haloes that is consis- 
tent with predictions of a standard ACDM cosmological model. The number density 
profile of the satellites appears to be shallower than of dark matter, with the scale ra- 
dius typically 1.6 times larger than of dark matter. The orbital anisotropy around red 
hosts exhibits a mild excess of radial motions, in agreement with the typical anisotropy 
profiles found in cosmological simulations, whereas blue galaxies are found to be con- 
sistent with an isotropic velocity distribution. Our new constraints on the halo masses 
of galaxies are used to provide analytic approximations of the halo-to-stcUar mass 
relation for red and blue galaxies. 
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1 INTRODUCTION 



According to the standard ACDM cosmological model, the 
majority of the total energy density of the Universe is 
deposited in the form of dark energy and dark matter 
(|Komatsu et al.l l201ll : IRozo et al.l I2OIOI '). The former is a 
homogeneously distributed component responsible for the 
observed acceleration of the Universe expansion, whereas 
the latter is highly clumped, setting up a base for the 
growth of cosmic structures. Dark matter (hereafter, DM) 
assembles within quasi-spherical haloes that host cosmo- 
logical objects of all scale, from dwarf galaxies to clus- 
ters of galaxies. The properties of such haloes are one of 
the most fundamental prediction s of the current cosmolog- 
ical model. As first discovered bv lNavarro. Frenk. fc White! 
(| 19951 . 1 19961 ). and confirmed in many more re cent and much 
better-reso l ved cosmological si mulations (e.g. ISpringel et al] 
I2OO5I . I2OO8I : iKlvpin et aLll201ll ). a key feature of DM haloes 
is the universal shape of their (hereafter NFW) density pro- 
file, whose logarithmic slope varies from —1 in the centre to 
—3 at large radii, while the transition scale between these 
two slopes is correlated with the halo mass (jNavarro et al.l 



I1997I : iLudlow et alll201ll ). This property is the subject of 
various observational tests at all halo masses. 

Studying the properties of DM distribution in galactic 
haloes is a challenge. Most methods rely on tracers whose 
positions coincide with the stellar component. Therefore, 
they probe only the inner part of an underlying gravita- 
tional potential of DM halo on the scale of a few per cent 
of the virial radius. Common means to study the inner 
part of DM density profiles is the measurement of rota- 
tion curves for spiral galaxies (jSofue fc RubirJ 200 j). the 
linc-of-sight vel ocity dispersio n profiles of stars ( Bertin et al.l 
1994 : Cappc Uari et al.| I2OO6I ) or planetary nebulae (e.g. , 
Napolitano et al. 201lh. strong (iKoopmans et af] 20061) 
and w eak lensing ( Mandelbaum et al.ll2006al: ICavazzi et all 
l200'iD or X-ray observations (|Humphrev et all |2006| ) for 
early type galaxies. The main difficulty in interpretation 
of these data arises from the fact that the mass of DM is 
comparable to the baryonic component and, therefore, con- 
straints on DM mass profile depe nds critically on the mas s 
estimate of the stellar component (|Mamon fc Loka^l2005al ). 
In particular, the uncertainty of stellar population models, 
when applied to elliptical galaxies, leads unavoidably to an 
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ambi guity about the shape of DM density profile (e.g.. lGrillol 
|2012| ). 

There are only two methods that allow to measure the 
DM distribution at distances comparable to the virial radius 
of the halo: weak lensing and kinematics of satellite galax- 
iesQ Due to a very weak signal per galaxy, both methods 
rely on stacking the data, giving insight into the properties 
of a spherically averaged rather than individual haloes. 

Lensing analyses were successfully applied to measure 
projected mass profiles in galacti c haloes. Most results are 
consi stent with a universal NFW (|Navarro. Frenk." fc White! 

density profile of dark matter and the mass- 
concentration relation emerging from cosmological simu- 
lations of a standard ACDM model (|Mandelbaum et al.l 
l2006al. I2OO8I'). Nev ertheless, one weak lensing study 
i Gavazzi etall |2007| ) concludes to a shallower DM den- 



sity profile at large radii around massive elliptical galaxies, 
even slightly shallower than the singular isothermal sphere 
model with p oc . Moreover, strong lensing studies of 
ellipticals also point to a density profile with slope very 
close to —2 between 0.3 and 0.9 Re. (jKoopmans et al"]|2006l . 
|2009| ')FI Howe ver, the naive superposition of NFW DM and 
the observed ISersid (| 19681 ) model for th e stars, leads to 
a slop e close to —2 (from the model of iMamon fc Loka^ 
l2005bl. we predict a slope of —2.2 ±0.1 in the range studied 
by iKoopmans et al.h: the stars domina te the mass profile 
within ~ 2Re { Mamon fc Lokasll2005bl ), but at large radii 
the NFW DM component should dominate and the slope 
should be considerably steeper than —2 (~ —2.6 at the virial 
radius from the model above). 

Satellite kinematics provide a popular means to esti- 
mate halo masses. Because host galaxies possess a very 
small number of observable satellites, usually one or two, 
one must stack the sa tellites over many host galax- 
ies. Still, early attempts JZaritskv. Smith. Frenk. fc White! 
! 19931 : !Zaritskv fc Whitell994 l suffered from their small sam- 
ple sizes. The Sloan Digital Sky Survey (SDSS) is the 
first very large spectroscopic sample of galaxies with accu- 
rate redshifts and digital p hotometry for a cre dible anal- 
ysis of satellite kinematics. !McKav et al.l (|2002l ) estimated 
host galaxy masses out to a fixed radius using M(rap) = 
Criipa^p{rs.p)/G, where C — dlnz^/dlnr is a constant de- 
termined by fitting a power-law to the stacked s atellite sur- 
face density profile. !Brainerd fc Specian! (|2003! ) perform a 
similar analysis on the 2dFGRS, where they measure the 
line-of-sight (LOS) velocity dispersion by fitting the LOS 
velocity distribution by a Gaussian plus a constant term for 
interlopers (instead of simply removing the high-velocity in- 
terlopers). They were the first to obtain M/L as a function of 
luminosity, and both for red and blue hosts. But their anal- 
ysis suffered from the relati vely inaccura t e velo cities and 
photometry of the 2dFGRS. !Prada et al] (|2003l ') were the 
first to notice a decline of LOS velocity dispersion at pro- 
jected radii. They showed that the distribution of satellites 
in projected phase spac e (PPS) is consistent with the ex- 
pectations from ACDM. !Conrov et al.! |2003) analysed the 



^ Strong lensing is typically restricted to inner regions, while X- 
ray measurements extend only to about half the virial radius. 
^ We call Re the effective radius, containing half of the luminosity 
in projection. 



satellites from the Data Release 4 (DR4) of the SDSS and 
from the DEEP2 survey ai z ~ 1, again with a model for the 
interlopers, and were the first to derive the variation of virial 
mass with host galaxy luminosity, separating red and blue 
galaxies. They found that red host galaxies of given blue 
luminosi ty have double th e halo mass as their blue coun- 
terparts. !More et al.l (|201l! ) added a second Gaussian to the 
Gaussian-|-flat distribution of LOS velocities and fit aper- 
ture velocity dispersions (which are less sensitive than LOS 
velocity dispersions to the unknown orbital anisotropy and 
its radial variation) to find a halo vers us stellar mass in close 
agreement with that of !Conrov et al.! . 

Unfortunately, all these analyses have flaws. For exam- 
ple, they all assume that the LOS velocity distribution is 
a Gaussian (generally plus a uniform interloper distribu- 
tion), while it is known that anisotropic ve locities lead to 
non-Gaussian LOS velocities (|Merrittl!l987^ . Thus by tak- 
ing into account the non- Gaussian nature of the L OS veloc- 
ity distribution (see also !Amorisco fc Evans! [20]! ). one can 
both obtain more accurate constraints on the mass profile 
and deri ve constraints on th e orbital anisotropy. 

We (!Woitak et al.|[2009! ) have recently developed a self- 
consistent method to derive at the same time the mass and 
velocity anisotropy profiles of spherical systems. Our method 
is based on the fac t that the distribution of o bjects in PPS 
is a triple integral (|Deionghe fc Merritt![l993 ) over the LOS 
and the two plane-of-sky velocities of the six-dimensional 
distribution function (DF) p arameterized in term s of energy 
and angular momentum that !Woitak et al.! (|2008h measured 
on the halos of a ACDM simulation. This approach gives 
much deeper insight into the data than tests of consistency 
shown before. It allows for a self-consistent comparison be- 
tween a set of physical parameters determined from cosmo- 
logical simulations and observations. Furthermore, analysis 
based on a PPS model does not rely on data binning which 
always introduces an artificial signal smoothing. 

This manuscript is organised as follows. In section 2, we 
describe the data and criteria for selecting isolated galaxies 
and their satellites. Section 3 presents our dynamical model 
and a method of constraining parameters of the systems. 
The results of data analysis and discussion are presented in 
section 4. The summary and a discussion follow in section 5. 
In this work, we adopted a flat ACDM cosmological model 
with Qm = 0.3 and Ho = 70 km s"^ Mpc"^ 



2 DATA 

We made use of th e Sloan Digital Sky Su rvey Data Release 
Seven (SDSS DR7. !Abazaiian et al.l!2009! ') to select isolated 
galaxies and the satellite galaxies orbiting them. To search 
for the host galaxies, we considered a volume-limited sub- 
sample of the spectroscopic part of the survey deflned by an 
r-band Petrosian absolute magnitude threshold Mr — —19.0 
and redshift range of 3000 km s"^ < cz < 25000 km s"\ 
The apparent magnitudes were converted to the absolute 
scale assuming fc-corr ection based on an analytical approx- 
imatio n provided by !Chilingarian. Melchior. fc ZolotukhinI 

iiolO). 

We defined isolated central galaxies as those that are 
brighter by AM then every other galaxy lying inside an ob- 
servational cylinder of a projected radius AR and a line-of- 
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sight velocity range 2Awios. We fixed aU fiducial parameters 
at values defining a rather restrictive criterion for galaxy 
isolation: AM > 1.505 (corresponding to the flux ratio 
of at least 2), AR < 1 Mpc and A^i o s < 1500 km 
(for comparison, see iMcKav et al.ll2002l: JPrada et al.ll2003l: 



■ den Bosch et al.ll2Q04l : IConrov et al.ll2007l : iKlvpin et al.l 
201ll ). All galaxies lying in the cylinder and that are dim- 
mer than the magnitude threshold are considered to be the 
satellites of the central galaxies. Due to a rather wide ve- 
locity cut-off, some of them are galaxies of background or 
foreground (interlopers). Disentangling between these two 
classes of galaxies is an intrinsic part of data analysis de- 
scribed in the following section. 

We split the sample of the selected host galaxies into red 
and blue galaxies usin g g — r + O.OlTMr colour diagnostic 
fsee lRoche et al.ll2010l ). where g and r are fc-corrected Pet- 
rosian magnitudes. A boundary value of this diagnostic was 
flxed at 0.25 which is a minimum of the colour distribution 
lying between two Gaussian components corresponding to 
two galaxy populations. Using the publicly available catalog 
of the stellar mass estimates from the SDSS DR7lf|we found 
the masses of the stellar component of all central galax- 
ies. Stellar m asses were estimat e d usin g Bayesian approach 
as outlined in lKauffmann et al.l (|2003|). with fitting th e ob- 
served photometry as described in lSalim e t al.' (200j). The 
model assumed initial mass function oflChabricr {20m). We 
neglected 8 per cent of host galaxies for which stellar mass 
estimates were not available. 

Our final sample consists of 8800 and 2600 satellites 
around 3800 red and 1600 blue hosts, respectively. The host 
galaxies cover the stellar mass range from logjg(M*/M0) = 
10.0 to logj^g(Af*/M0) = 11.8 for red galaxies and from 
logio(M^/Mo) = 9.5 to logio(M^/Mo) = 11.0 for the blue 
ones. Since the stellar mass is a better indicator of the haloes 
mass than galaxy luminosity (|More et al.ll20lil ). we split the 
sample of the host galaxies into several bins of the stellar 
mass. We used 6 and 3 bins for the red and blue hosts, 
respectively, as indicated in Table [T] This procedure guar- 
antees that the kinematic sample in every bin represents a 
homogeneous sample of DM haloes. 



3 INFERENCE OF PHYSICAL PARAMETERS 

The velocity distribution of satellite galaxies is mostly 
determined by the gravitational potential of DM haloes 
of the central galaxy. The second factor is the orbital 
anisotropy describing the fraction of radial-to-tangential or- 
bits in the system. This additional degree of freedom makes 
data analysis more complex due to a well-known fact of 
the mass-anisotropy de generacy (|Binnev fc MamonI 1 19821 : 
iMerrifield fc Kent|[l990l ) . Breaking this de generacy requires 
using rather complicated models accounti ng for higher-order 



corre c tions t o the J eans equation (e.g.. IMerrifield fc Kent 
1990l : iLokad |2002| : iLokas fc MamonI l2003l : IWoitak et al 



20091 ). On the other hand, the advantage is that the same 



data allow to study two physical properties of the host- 
satellites systems at the same time - mass distribution of 



DM halo and the orbital structure of the satellites (e.g . , 
iLokas fc Mamonll2003l : lLokasll2009l : IWoitak fc Lokasll2010l ). 

We analysed the kinematic data in terms of the pro- 
jected phase-space density, i.e. the density of satellite galax- 
ies on the plane spanned by the line-of-sight velocity vios 
and the projected distance from the host galaxy R. Due to 
projection effects, the vios — R plane is populated by the true 
physical satellite galaxies of the central galaxies as well as 
interlopers. We did not apply any interloper removal to the 
data, but rather we accounted for the presence of interlopers 
in a statistical sense as an inherent part of a proper analysis. 
This approach is particularly justified in case of the compos- 
ite kinematic data for which the phase-space distribution 
of interlopers is smooth enough to be modelled by a con- 
tinuous probability function, and is commonl y adopted in 
many studies on kinemat i cs of satellite galaxies JPrada et al.l 
I2OO3I : IConrov et aflbOOTl : iKlvpin fc PradalbOOgf y The prob- 
ability describing the phase-space distribution of interlopers 
introduces an additional degree of freedom to the proper 
model associated with the physical properties of the host- 
satellites systems. The observed projected phase-space den- 
sity pios(-R, wios) may be expressed as the following sum 



g{R,Vlos) = (1 - Pi) gs{R,Vloa) +Pi 5i(-R, t^los) : 



(1) 



where gaiRjVios) and (/i(ii, uios) are the projected phase- 
space densities of satellite galaxies and interlopers respec- 
tively, and Pi is the probability of a randomly picked galaxy 
being an interloper. 

3.1 Phase-space density model 

As a model of the phase-space distribution of satellite galax- 
ies, we used an anisotropic model of the distribution func- 
tion developed by (jWojtak ct al. 2008). The model was de- 
signed to describe the phase-space properties of simulated 
DM haloes and it was successfully utilised to constrain mass 

J irofile and the orbital anisotropy in nearby galaxy clusters 
Woitak fc Lokasll2010l ). The only modification required to 
adjust the model to the new context of systems of hosts 
and satellites is a non-constant ratio of DM-to-tracer den- 
sity profile. Constant mass-to-light ratio appears to be a ro- 
bust assumption in gal axy clusters (|Biviano fc Girardi|[200l : 
iLokas fc MamonI |2003| ). but it is not justified for the sys- 
tems of satellite galaxies for which observations point to a 
bias between sp atial distribution of DM and the satellite 
iGuo et al.ll201S 



Following IWoitak et aU 120081 ). we considered the 
phase-space density /(r, v) of the following form 



/(r,v) = /B(i5)L 



-2/3o 



2Ll 



(/3o-/3oc) 



(2) 



where E and L are positively defined binding energy and an- 
gular momentum per unit mass, /?o and Poo are the asymp- 
totic values of the anisotropy parameter at small and large 
radii, respectively. The anisotropy parameter quantifies the 
orbital anisotropy in terms of the ratio of the radial-to- 
tangential velocity dispersion and is traditionally defined as 



/3(r) 



1 



2 0-2 



(3) 



where a^- and at are the velocity dispersions in radial and 
http://www.mpa-garching.mpg.do/SDSS/DR7/Data/steIIarmass.htrrtiangential direction, respectively. The angular momentum 
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part of the distribution function ((2)1 is a generalisation 
of a well-know L~'^^ ansatz for systems with a constant 
anisotropy parameter (|Hcnon 19731 ) . It permits a wide family 
of the anisotropy profiles, suitable for constraining not only 
a global degree of the anisotropy, but also its radial profile. 
The anisotropy profiles are monotonic functions changing 
between two asymptotic val ues with the radius of transition 
given by Lo parameter (see lWoitak et allbOOSl . for details). 

The energy part of the distribution function ((2)1 is re- 
lated to the density profile of satellite galaxies ps (r) and an 
underlying absolute value of the gravitational potential ^'(r) 
through the following integral equation 



/s(i5)i-^'^"(l + A^ 



(,8o-/3oc) 



d^«. 



(4) 



where E = ^{r) — ^v^ is the binding energy. The mass profile 
at large distances from the central galaxies is dominated by 
DM. While the stars can contribute up to half the mass 
at 5 R e (see upper right panel of Fig. 4 of iMamon &i Lokaj 
l2005bl ). we choose to neglect the contribution of stars (or 
in other words to assume that they are distributed like the 
DM). We checked that this assumption has negligible effect 
on our analysis (see Discussion). 

We also approximated th e DM density profi le by 
the universal NFW profile (jNavarro et al.1 Il997l ) for 
which the gra vitation al potentia l takes the following form 
l|Cole fc Lacev 1996; Lokas fc Mamon..2001i ') 



*(r) = *o 



ln(l + r/rP^) 



(5) 



where is the scale radius at which logarithmic slope of 
DM density profile equals to —2. Our choice of the NFW pa- 
rameterisation is motivated not only by cosmological simula- 
tions, but also by observational results showing consistency 
between satellite kinematics and dynami cal predictions for 
DM haloes with th e NFW density profile (|Prada et al.ll2003l : 
iKlvpin et al]|201ll '). We note that the *o and are the 
principal parameters in the analysis which can be easily con- 
verted into more popular quantities describing the mass pro- 
file of DM haloes s uch as the virial mass o r the concentra- 
tion parameter fsee lLokas fc MamonlliooH for all equations 
needed for parameter transformations). 

The number density profile of satellite galaxies may be 
effectively approximated by the NFW profile with a scale ra- 
dius unrelated to the concentration of DM (e.g., IGuo et al] 
|2012| ). This property is independent of morphological type 
and the stellar mass of the host galaxy. Following this ob- 
servational motivation, we adopted 



Ps(r) oc 



(6) 



as the number density of the satellites, where Tg'** is a new 
scale radius. 

Having specified ps{r) and ^{r) one can solve equa- 
tion Q fo r the energy p a rt of t he distribution function. As 
shown by IWoitak et al.l (|2008[ ). the integral over velocity 
space may be reduced to a one-dimensional problem. Then, 
the resulting integral equation may be inverted numerically. 
We use d the same scheme o f the integral inversion as out- 
lined in IWoitak et al.l (|2008l ) . Although the algorithm was 
designed to work for a single-component system with an 



NFW density profile, we checked that it is also feasible in 
case of two-component systems. 

The projected phase-space density of satellite galax- 
ies in ^ was obtained by integrating the full phase-space 
density ([2} over velocities vj_ perpendiculax to the line of 
sight and a spatial coordin ate z parallel to the line of sight 
(|Deionghe fc Merrittiri992l) 



g{vios,R) = 2tvR / dz 



J{E,L). 



(7) 



Following the scheme outlined bv lWoitak et al.l (|2009l ). we 
calculated this integral numerically using Gaussian quadra- 
ture. We did not apply any fiducial truncation to the dis- 
tance along the line of sight keeping the upper limit of the 
corresponding integral as defined by the condition of positive 
binding energy, i.e. E = *(V-R^ + z^) - I'L/S ^ 0. 

For interlopers, we adopted a uniform projected phase- 
space distribution, i.e. oc R. It has been demonstrated 
that this model effectively separates gravitationally bound 
members of the central system from interlopers whose veloci- 
ties are mostly dominated by the Hubble flow jWoitak et al.l 
l2007l : iMamon. Biviano. fc Murantell20ldl ) . 

As a consistency check, we verified the robustness of this 
model by comparing it with the distribution of galaxies at 
R > 100 kpc and 1000km s"^ < \vias\ < 1500kms"\ which 
constitutes a clear subsample of interlopers for all bins of 
the stellar mass (|vioh| exceeds maximum escape velocity at 
all halo masses as large as lO^^'^M©). Although the surface 
density of these galaxies is consistent with a uniform distri- 
bution, the velocity distribution for massive red hosts ap- 
pears to exhibit a weak and broad non-uniform component. 
We attribute this feature to an environmental effect: massive 
red galaxies, although selected by the same isolation crite- 
ria, tend to populate denser environment which enhances 
the number of interlopers relative to low-density environ- 
ments behind and in front of the host galaxy. Consequently, 
interloper density decreases with increasing velocity \vioe\- 
We found that in order to account for this effect it suffices 
to consider the following modification of a uniform model 



giocR \ (1 -pg) +pg 



2'Kaa 



■ exp 



2^1 



(8) 



where erg = 500 km s~^ and pg is a nuisance parameter 
defining relative weight of the Gaussian component. This 
form of gi was used only for three bins of the most massive 
red host galaxies with logj^Q(M*/MQ) > 11.0. 

3.2 Incompleteness 

Spectroscopic survey of the SDSS is not complete on angular 
scales smaller than 55" imposed b y the minimum sepa ration 
of the fibers in the spectrograph (jBlanton et al hoO^ ). The 
limit of the completeness corresponds to the physical scale 
of 92 kpc at the maximum redshift defining the sample of 
isolated galaxies. This distance is a substantial fraction of 
the virial radius and, therefore, it is relevant to correct the 
projected phase-space density model for the incompleteness. 
The correction was incorporated by means of weighting the 
projected phase-space density according to the local value 
of the completeness 'w{R) 



g{R,vios) = [g {R,vios)] 



l/m(R) 



(9) 
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10.0<log Mj,<11.0, red 
11.0<log Mj,<11.8, red 
9.5<log M„<10.5, blue 
10.5<log M„<11.0, blue 



R [kpc] 

Figure 1. Completeness of the SDSS redshift survey around iso- 
lated host galaxies. The data points are the measured values of 
completeness, obtained by comparing the counts of SDSS spec- 
troscopic and photometric data. The colours represent different 
stellar mass bins of red and blues host galaxies as indicated in 
the legend. 



We measured the completeness of the data by comput- 
ing the ratio of the surface number density of galaxies with 
spectroscopic redshifts to the surface density of all galaxies 
brighter than the magnitude limit of the SDSS spectroscopic 
survey, i.e. r = 17.7 in r-band of Petrosian magnitude (see 
Fig. [!}. In order to account for a redshift dependence we 
split the sample of red and blue galaxies into two bins of 
the stellar mass corresponding to two classes of the absolute 
magnitude. We found that the completeness may be well 
fitted by the following analytical form 



w{R) 



: WO- 



[R/ Rv. 



(10) 



1 + {R/R^)°' ' 

where wq, R„ and a are free parameters. We fitted this pro- 
file to the data in all mass bins of the host galaxies. The 
resulting best-fit profiles of w{R) (see solid lines in Fig. m 
were used as the final weighting function in ®. The final 
procedure of parameter inference was positively tested on 
incomplete mock data generated from the distribution func- 
tion with a uniform background of interlopers. 

3.3 Parameter estimation 

We made use of kinematic data of satellite galaxies to place 
constraints on parameters of DM mass profile and the or- 
bital anisotropy. For this, we adopted a Bayesian approach, 
maximising the likelihood of the distribution of satellites in 
PPS. The likelihood function £ was defined as 



ln£ = Y.^ng{Rj 



E 



-w{R, 



■ In g{R, 



!a) 



(11) 



where the sum is over all satellite galaxies in a given stellar 
mass bin of the host galaxies and a is a vector of the model 
parameters. Both gs and gi are normalised to 1 over the 
area confined by velocity cut-off |uios| < 1500 km s~^ and 
radius range [-Rmin , i^max] . As the maximum radius Rmax, 



we used the virial radius (see next paragraph) as it de- 
termines a natural boundary of the equilibrated part of 
DM haloes. Since the virial radius is a function of some 
model parameters, its value was estimated in an iterative 
approach starting with the best initi al guess based on th e 
halo-stellar- mass relation provided bv lDutton et al.l (|2010l ). 
For the most massive host galaxies, we imposed an addi- 
tional limit i?max < 400 kpc which prevented from includ- 
ing the satellites which may be common to the host galaxy 
and local group or cluster of galaxies. We also imposed a 
minimum radius, because 1) our correction to the spectro- 
scopic incompleteness is uncertain at small radii where our 
measured completeness is low, and 2) because photometric 
pipelines such as that of the SDSS tend to fragment large 
(host) galaxies into one big one and many small ones sur- 
rounding it, that appear like satellites, but are HII regions or 
spiral arms instead. The minimum radius -Rmin was fixed at 
5 effective radii for red galaxies and 15 kpc for the blue ones. 
Effective radii were estimated independently for every stellar 
m ass bin using a scal i ng rel ation with the stellar mass found 
bv lHvde fc Bernardil (|2009l '). The resuhing minimum radius 
Rmin changes from 10 kpc for logio(M^/Mo) = 10.0 - 10.5 
to 55 kpc for logio(M*/Mo) = 11.5 - 11.8. 

The ma:ximisation of the likelihood was carried out us- 
ing the Markov Chain Monte Carlo (MCMC) technique w ith 
the Metropolis-Hastings algorithm (|Gelman et al.|[2004l . see 
e.g.). The set of the primary parameters used in the MCMC 
analysis comprises the scale radius Tg"* of the number den- 
sity of satellites ((6|, the dimensionless r^^^/r^'^ ratio, the 
normalisation "to of the gravitational potential the 
asymptotic velocity anisotropy at small and large radii (/3o 
and Poa, respectively), the interloper probability pi in ^ 
and the relative weight pg of the Gaussian part in the veloc- 
ity distribution ((S]) of interlopers (applied only to the data 
of red galaxies with > lO^^M©). We determined con- 
straints on the mass profiles by converting parameters of 
the gravitational potential, i.e. ^'o and r^^, into the stan- 
dard parameters characterising DM halo with the universal 
NFW density profile; the virial mass Ma and the concentra- 
tion parameter ca- The virial mass is defined in terms of the 
mean density inside the sphere of radius ta (the so-called 
virial radius) relative to the critical density pc 



3Ma 
47rr^ 



(12) 



where A is the virial overdensity. The concentration param- 
eter is the virial radius expressed in the unit of the scale 
radius r^^, i.e. ca = rA/f?^. We adopted two commonly 
used values of the overdensity parameter: A = 200 and 
A — 100. The latter corresponds, to a 3 per cent precision, 
to the virial overdensity of a st andard ACDM cosmological 
model iBrvan fc NormanlllQOsI '). 

We carried out the MCMC analysis assuming log- 
uniform priors for r^^', and r^^^ /r^^, and uniform pri- 
ors for all remaining parameters. We fixed Lo (eq. [2]) at 
0.2v'Wors^' corresponding to the ~ 1 r^'^'^ transition radius 
bet ween two asymptotic v alues of the anisotropy parame- 
ter (|Woitak fc Lokail2010l V The central anisotropy po was 
limited by /3o < 1/2 in order to the prevent distr i butio n 
function from taking negative values (|An fc Evans! |2006| ) . 
Constraints on all parameters are based on Markov chains 
containing 2 x 10* models in every bin of the stellar mass. 
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Table 1. Constraints on the parameters of DM density profile and the number density of galaxy satellites in all bins of the stellar mass 
and both types of host galaxies (red and blue): DM halo mass Mioo (or M2oo)i the concentration parameter cioo (or C200) a-nd the ratio 
of traccr-to-DM scale radius ''g(sat)/''s(DM) ' "^^^ table provides the best fit values at the maximum of the posterior probability and 
the ranges containing 68 per cent of the corresponding marginal probability. The last column contains the number of satellites between 
Rniin and Rmax, A^sat , and the number of hosts, N-^^ost- 
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1 I ^^^^ ' ' — ' ' — ■— ^ 

0.8 1 2 4 6 8 10 20 40 60 

MiooLIO'^Mq] 



Figure 2. Concentration-mass relation for galaxy-size haloes. 
The points with the error bars show constraints on the halo mass 
and the concentration parameter of the DM density profile in- 
ferred from the kinematics of the satellite galaxies. The purple 
line is the best-fit power-law fit and the shaded region is the 
1(T confidence area resulting from propagating the errors on the 
slope and the normalisation (dotted lines show the best-fit profile 
and the la confidence area when neglecting the most outlying 
point at Afioo = 1-5 X lOl-^M©). The black solid line is the me- 
dian concentration-mass relation from cosmological simulations of 
a sta ndard ACDM model (the Bolshoi Simulation, iKlvpin et alj 
[20O). 



Every chain was preceded by a number of trial chains ran to 
estimate the covariance matrix o f the proposal probability 
distribution (jCelman et al.ll2004l '). 



4 RESULTS 
4.1 Mass profile 

Table [1] contains constraints on the halo mass, the concen- 
tration parameter and the ratio of the tracer-to-dark-matter 
scale radius of the density profile, for all bins of the host 
galaxies. Independently of the halo mass and the colour of 
the hosts, the concentration parameter appears to be fairly 
well constrained with a typical value ~ 10. This is a plausi- 
ble signature that DM density profile is steeper then r^^ at 
large radii [r > rioo/10) and shallower at small radii. 

Fig.[2]shows the comparison between our 'observational' 
constraints on the concentration-mass relation from satellite 
kinematics (red and blue points) and its theoretical pre- 
dictions (black solid line). For the latter, we plotted the 
concentration-mass profile from the Bolshoi Simulation - a 
high resolution simulation of a standard ACDM cosmolog- 
ical model with the most updated cosmological parameters 
(|Klvpin et al.|[201ll ). The concentration parameters inferred 
from the satellite velocities are fairly consistent with the pro- 
file from cosmological simulations. However, our constraints 
are not tight enough to determine robustly the slope of the 
mass-concentration relation. Power-law fits to the data of 
both types of the host galaxies, with symmetrised errors, 
yields the slope —0.38 ± 0.22 (purple solid line in Fig. [2]) or 
—0.16 ± 0.22 (purple dashed line in Fig. [2| when neglecting 
the most outlying point at Mioo ~ 1-5 x lOi ^Mp. This is 
consi stent with the predicted value of —0.075 (|Klvpin et al.l 
I2OIII) as well as with a flat profile. The best-fit normalisa- 
tion at Mioo = 3 X lOi^M0 is 8.6 ± 1.8 when fitting to aU 
data points, and 9.1 ± 1.6 when excluding the outlier, in 
excellent agreement with the concentration Cioo =9.1 from 
the simulations of the current ACDM cosmological model 
iKlvpin et al.ll201ll ). 

Since the data appear not to reveal a statistically sig- 
nificant dependence of the concentration parameter on the 
halo mass, we combined constraints on the concentration 
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Figure 3. Probability distribution for the concentration parame- 
ter of DM density profile, as combined from all bins of the stellar 
mass. The arrows indica te the median conc entration parameters 
of simulated DM haloes jKlvpin et alJl2011^ corresponding to the 
median halo mass in red, blue and both types of the central galax- 
ies. Red, blue and black colours correspond to the red and blue 
host galaxies, and both types of the central galaxies. 



from all bins of the stellar mass. Fig. [3] shows the result- 
ing posterior probability distribution for red, blue and both 
types of the host galaxies. The concentration parameter of 
DM haloes hosting red galaxies is fully consistent with the 
predictions of DM cosmological simulations corresponding 
to the same range of the halo masses (red arrow). DM den- 
sity in haloes of blue galaxies appears to be slightly less 
concentrated compared to theoretical expectation (blue ar- 
rows). This deviation, however, is not statistically significant 
(around la) and it may be easily verified in future with the 
use of more data from upcoming redshift surveys. Combin- 
ing results from all bins of the stellar mass yields the typical 



concentration parameter 8.OI2 g for red galaxies, 5.7'!i2't for 
the blue ones and 6.8l^ g for both types. 

Fig. |4] shows probability distribution for rl'^^/r^^ com- 
bined from all bins of the stellar mass. Typical scale ra- 
dius of the satellite number density profiles around red host 
galaxies is larger by factor of l.Slp g than that correspond- 
ing to DM (r^'"/^?^ > 1-0 at the 97 per cent confidence 
level). This finding confirms previous studies based on pho- 
tometric data from the SDSS and showing that the spatial 
distribution of the satellite around red galaxies is more ex- 
tended than of DM and is well-fitted by the NFW profile 
with the concentration par ameter typically 2 times smaller 
than that expected for DM (|Guo et al.l2012l ) . Constraints on 
the ratio of the sateUite-to-dark-matter scale radius for blue 
galaxies do not point to any bias between the spatial distri- 
bution of DM and the satellites, i.e. r^Vr?^ = 1.4l°;*. 
Combining the results from both types of central galax- 
ies yields the host-independent ratio rrV*"?^ = l-Sto.s 
(rrV^'s^^ > 1.0 at the 97 per cent confidence level). This 
estimate is fully consistent with the bias between the con- 
centration parameters of the subhalo number density profile 
and dark matter density profile measured in cosmological 
simu lations, ff'^^/rP^ = 1.5 from the Millennium Simula- 

1.3 from the Bolshoi 



tion (ISales et al.ll2007^ arid r^V rP^ 
Simulation (|Klvpin etHlbOllI ') . 



Figure 4. Probability distribution for r^^^ /r^^ combined from 
all bins of the stellar mass, where r^^' and are the scale radii 
of the satellite number density profile and dark matter density 
profile. Red, blue and black colours correspond to the red and 
blue host galaxies, and both types of the central galaxies. 



4.2 Halo-stellar mass relation 

Fig. [5] shows the halo mass as a function of the stellar mass 
of both types of host galaxies. The halo mass is measured 
with accuracy up to 0.05 dex for red galaxies with stellar 
masses logjQ(M*/M0) = 10.75 — 11.25 for which number 
of the satellites per stellar mass bin reaches maximum. Ob- 
tained constraints on the halo-stellar mass relation reveal 
a change of the slope between low and high stellar masses, 
with the transition mass log]^Q(M*/M0) w 11. Following 
iDutton et al] (|2010l ). we find that a reasonable fit to the 
data is achieved using the following function 



M2, 



(b-a)/-, 



(13) 



where a and /3 are the logarithmic slopes at small and 
large stellar masses, respectively, M*o and Mho are the 
stellar and halo mass at the transition point, and 7 is a 
parameter controlling the sharpness of the transition. Fit- 
ting this function to the data of red galaxies yields a = 
0.29, b = 2.91, 7 = 1.24, logi(,(M*o/M0) = 11.3 and 
logj^Q(M/io/MQ) — 13.1. Constraints on the halo-to-stellar 
mass relation for blue galaxies cover approximately one or- 
der of magnitude in the stellar mass and are consistent with 
a power-law profile with the logarithmic slope 0.66 ± 0.07 
and the normalisation logj^Q(M2oo/MQ) = 12.0 at the stellar 
mass logio(M*/M0) = 10.3. 

We compare our re sults with empirical profiles obtained 
bv lDutton et al.l (|2010l ) as a compilation of all observational 
constraints on the halo-stellar mass relation available in the 
literature (see shaded stripes in Fig. [5]). The widths of the 
profiles correspond to the 2a credibility range, where a is 
the error dominated by systematics related to differences 
bet ween various mass m easurements based on weak lens- 
ing (|Mandelbaum et al.|[2006h. ,2008: Schulz ct al. 20101) o r 
satellite kinematics (|Conrov et al.ll2007l : [More et aLlbOllI ). 
and dashed lines represent the best fit profiles. We find that 
the halo-stellar mass relation from our analysis of red galax- 
ies is fairly consistent with other measurements. It only 
exhibits a sharper transition between low and high stellar 



© 0000 RAS, MNRAS 000, 000-000 



8 



R. Wojtak & G. A. Mamon 




10.5 11.0 
loglo(M**lo) 



12.0 



1.0 

0.8 
0.6 
0.4 
0.2 
0.0 
-0.2 
-0.4 
-0.6 
-0.8 
-1.0 















■ 

/\ , . <?: 





-1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 



Figure 5. Relation between the halo mass M200 and the stellar 
mass M* for rod and blue galaxies (red and blue colour, respec- 
tively). The error bars indicate the range containing 68 per cent 
of the marginal probability distribution for the halo mass and 
the scatter of the stellar masses in bins. Solid lines are the best 
fit power-law profile for blue galaxies and the best fit double- 
power-law profile given by I I13I I for red galaxies. Shaded stripes 
and dashed lines show halo-stellar mass relation r epresenting a 
comp ilation of results available in the literature jPutton et alj 
l201(t) . The widths of the stripes correspond to the la credibility 
range and the dashed line is the best fit double-power-law profile. 



mass regime and a slightly shallower profile at small stellar 
masses. 

The halo-stellar mass relation for blue galaxies appears 
to be offset by 0.15 dex towards higher halo masses, al- 
though the trend is not statistically signific ant. We note, 
howe ver, that our reference profile plotted (|Dutton et al.l 
represents a compilation of a quite heterogeneous 
sample of measurements which are incidentally obtained 
m ostly from satellite k inematics. Closer inspection of Fig. 1 
m button et aD (j201(t ) shows that most of the halo masses 
lie above the best-fit profile. In part icular, halo masse s ob- 
tained by IConrov et al] (|2007l ) and iMore et al.] (|201ll ') are 
offset by 0.20 dex and 0.1 dex respectively which are suffi- 
ciently large values to match our measurement of the halo 
masses. It appears t hat the halo-s t ellar mass relation for 
blue galaxies found bv lDutton et al.l ()2010l ) is biased towards 
lower halo masses, probably as a result of special selection 
of data points used for fitting. 



4.3 Anisotropy of the satellite orbits 

Constraints on the anisotropy parameter obtained in in- 
dividual bins of the stellar mass are not tight enough to 
draw any solid conclusion. For example, working in separate 
bins of host stellar mass, we cannot differentiate between 
an isotropic velocity distribution (/3 = 0) and the typical 
anisotropy profile found in simulated DM haloes where /3 in- 
creases with radi us from 0.1 i n the halo centre to 0.3 — 0.5 at 
the v i rial radius ( Woitak et al..,2005. : .Ascasibar fc Gottlobeij 
l2008l : ICuesta et al.ll20o5 r Since there is no theoretical hint 
that the anisotropy may depend on the halo mass and re- 
sults obtained in different stellar mass bins do not reveal 
any trend with the stellar mass, it is advisable to combine 
constraints from all bins into one. 



Figure 6. Constraints on the asymptotic values of the anisotropy 



profile parameter /3(r), /3o = {^{t ^ r^'^*) and fi^ 



{r > r-|^'). 



The contours are the confidence regions containing 68 and 95 per 
cent of the probability combining results from the stellar mass 
bins. Red, blue and black colours correspond to red, blue and 
both types of host galaxies. 



Fig. |6] shows contours of the resulting joint probability 
distribution of two parameters determining the anisotropy 
profile, Po = P{r < rf*) and /3oo = /3(r > rf*). When we 
combine all host galaxy mass bins, we find that the satellite 
orbits are mildly radially anisotropic with Pq — 0.2 ± 0.1 
and /3oo = 0.3 ± 0.2, in fair agreement with the predictions 
from cosmological simulations as well as with the universal 
relations between the anisotropy an d the logarithmic slope 
of an underlying DM density profile ([Hansen fc Moordliooel : 
[Hansen et al-llioiol ). The shift between the marginal proba- 
bility distributions for /3oo and Pq may indicate a tendency 
for the anisotropy to increase with radius. This effect, how- 
ever, is of a marginal statistical significance and the data 
still permit a flat anisotropy proflle. 

Our measurement of the anisotropy parameter relies 
mostly on the data for red galaxies (compare the red and 
black probability distributions in Fig. (6)1 . Combining results 
for red galaxies yields fio = 0.26lE!;?^ and = O-SStolr 
which are statistically consistent with the constraints ob- 
tained for both types of central galaxies. Measurement of 
the orbital anisotropy around blue galaxies is much weaker 
due to significantly smaller number of the satellites and it 
does not allow to differentiate between an isotropic velocity 
distribution and the anisotropy profile motivated by cosmo- 
logical simulations (see the blue probability distribution in 
Fig-O- The analysis of the combined probability distribu- 



tion results in /3o 



-0 2+° * 
'-'■^-0.6 



and /3a 



0.2 



-1-0.3 



5 SUMMARY AND DISCUSSION 

We made use of satellite kinematics compiled from the red- 
shift catalog of the SDSS to measure the physical proper- 
ties of dark matter haloes and the orbital anisotropy of the 
satellites orbiting fairly isolated galaxies. Data analysis was 
carried out in the framework of the projected phase-space 
density based on an anisotropic model of the distribution 
function for equilibrated spherical systems. This approach 
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avoids arbitrary binning the data. It furtiiermore allows us 
to break the mass-anisotropy degeneracy. 

The concentration parameter of DM density profiles has 
a typical value of ~ 9, in full consistency with th e results 
from cosmological simulations (|Klvpin et al.llioill ). We did 
not find statistically significant differences between DM con- 
centrations in red and blue galaxies. Our constraints on the 
concentration-mass relation are not tight enough to deter- 
mine its slope over the range of galactic halo masses. How- 
ever, a more robust measurement of the c — Mioo slope may 
be achieved by means of complementing these results by 
similar constraints at higher halo masses available in the 
literature. Combining the normalisation of the c — Mioo 
resulting from s a tellite kinematics with that obtained by 
IWoitak fc Lokaj (|2010l ) from galaxy kinematics in clusters 
(cioo = 6.8 ± 0.7 at Mioo = 4.9 x IO^Mq) yields the slope 
—0.05 ±0.04, in fair agreement with —0.075 from cosmolog- 
ical simulations (|Klvpin et al.|[201ll ). 

Although our analysis relies on a specific parameterisa- 
tion of the density profile, the existence of a well-constrained 
dark matter scale radius < rioo suffices to conclude 

that the observed satellite kinematics is fully compatible 
with the NFW density profile of DM and can hardly be rec- 
onciled with an isothermal sphere model suggested in sev- 
eral stud ies based on lensing analysi s of massive elliptica l 
galaxies (|Koopmans et al" l l2006l . I2OO9I : iGavazzi et aLll2007l ). 
An isothermal density profile of DM would noticeably affect 
the measurement of the concentration parameter resulting 
in « rioo. This finding confirms and complements a 

number of tests showing consistency of satellite kinemat- 
ics with the NFW pro file of DM density (|Prada et al.ll2003l : 
iKlvpin fc Pradall2009^ . It also contradicts results of some 
studies based on weak lensing analysis around massive el- 
liptical galaxies and resulting in isothermal sphere density 
profile at radii dominated by DM, i.e. ~ 100 effective radii 
(see e.g. iGavazzi et al.|[2007t ). This apparent flattening of 
the density proflle should probably be attributed to a pro- 
jection effect of the local dense environment rather than to 
DM haloes. 

In our analysis, by adopting a single NFW model for the 
host mass distribution, we have neglected the contribution 
of the stellar component. Since elliptical galaxies are known 
to be domi nated by their ste llar com ponent within the effec - 
tive radius (jMamon fc Lokaai,2005a : .Humphrev et al.ll200d ) , 
one may worry that our DM concentrations will be overes- 
timated, even though we only considered satellites further 
than 5 Re from the host galaxy. For example, the total den- 
sity proflle of a two-component NFW-(-Sersic model is very 
close to a singular isothermal in the range 0.1 — » 1 Re (see 
upper left panel of Fig. 4 of iMamon fc Lokad l2005bl ). as 
also determined from st rong lensing analyses i n a similar 
range of projected radii jKoopmans et al]l2006l ). Since the 
mass distribution returned from our model is most sensitive 
to the line-of-sight velocity dispersion profile, aiosiR), we 
asked ourselves how much lower would the DM concentra- 
tion be if we incorporated a lSersid (| 19681 ) model to the mass 
distribution. 

We performed this test for our red host galaxies. For 
each of our bins of stellar mass, we determined the median 
r-band absolute ma gnitude, then obtained from Table 3 of 
ISimard et al.l (|201lh the effective radii. Re, and Sersic in- 
dices, n, for these absolute magnitudes, after restricting the 



SDSS sample of ISimard et al.l to the Red Sequence (using 
the same cut as we did in our mass analysis) and our adopted 
redshift range. We then computed trios (ii) for the satellite 
population with a two-component NFW-I-Sersic mass model, 
adopting the central stellar mass of our mass bin with the 
values of Re and n that we obtained above, the satellite 
scale radius that we previously measured (derived from Ta- 
ble[T]), adopting /3(r) = 0.2 for the satellites (consistent with 
Fig. 6), as well as the DM normalisation derived from our 
1-component mass model (Table [l|. The DM concentration, 
cioo, is a free parameter. We iterated on cioo until our profile 
of aios{R) matched the profile expected for a 1-component 
NFW model. For each of the six mass bins, we were then 
able to match aic,s{R) to better than 0.5 per cent between 
5 -Re and the virial radius rioo. In the end, we found that 
the concentration of the DM component was 7 to 20 per 
cent lower than in the 1-component model. This suggests 
that our derived values of cioo are overestimated by 7 to 20 
per c ent. Note that if we allowed for adiabatic contraction 
(e.g., iGnedin et al.llioill ). the DM concentration would be 
greater, hence our overestimate would be lower. This sim- 
ple analysis suggests that our choice of a single NFW model 
for the mass distribution of isolated galaxies beyond 5 Re is 
reasonable, as the predicted line-of-sight velocity dispersion 
profiles match that of more realistic two component models 
and our concentration is typically overestimated by only 10 
per cent. 

Satellite kinematics reveals a bias between the spatial 
distributions of DM and satellite galaxies. The scale radius 
of the satellite number density profile is typically larger by 
factor of 1.6 than of DM. This finding is consistent with 
the estimate of a counterpart bias between DM particles 
and subhaloes found in cosmologica l simulations (see e.g. 
ISales et al.ll2007l : iKlvpin et all 1201 ll ). We did not find any 
statistically significant difference between the bias of red and 
blue galaxies. 

The orbital anisotropy of satellite galaxies exhibits a 
mild excess of radial orbits with typical anisotropy parame- 
ter P — 0.2 ±0.1 in the inner regions and /? — 0.3 ± 0.2 in 
the outer regions of their hosts. These constraints on the in- 
ner and outer asymptotic values of the anisotropy profile are 
too weak to reveal any statistically significant trend of the 
anisotropy with radius. Due to the substantial difference be- 
tween the numbers of satellites around red and blues hosts, 
our constraints on the orbital anisotropy come principally 
from the satellites around red hosts. This radial anisotropy 
around red (giant elliptical) galaxies has been pr edicted from 
hydr odynamical simulations of binary mergers (|Dekel et al.l 
I2OO5I ). Some giant ellipt ical galaxies shows signs of such 
radial outer anisotropy (iDas et al.l 20081: de Lorenzi et al] 
l2008t ). while other do not (|Napolitano et al.l 1201 ll ). In fact, 
our asymptotic value for the anisotropy. Poo = 0.3±0.2, sug- 
gests that satellites have less radial orbits than the particles 
at analogous radii within ACDM haloes. The velocity dis- 
tribution of the satellites orbiting blue galaxies is consistent 
with an isotropic model. 
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